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SYNOPSIS 


.V 


The  objective  of  this  program  is  to  experimentally  examine  the 
interaction  of  a  turbulent  mean  flow  with  entrained  air  bubbles. 
Particular  attention  has  been  paid  to  the  determination  of  the  relative 
acoustic  radiation  efficiencies  of  bubbles  undergoing  simple  harmonic 
oscillation  and  those  undergoing  splitting.  There  were  three  phases  to 
this  program: 

1)  Design  and  fabrication  of  a  bubble  Injection  system. 
Redesign  and  quieting  of  a  turbulent  jet  facility. 

2)  Design  and  fabrication  of  a  facility  to  study  bubble  for¬ 
mation  noise.  Acoustic  and  high  speed  observations  were  made 
and  correlated  with  new  theory  developed  in  this  program. 

3)  Experimental  determination  of  the  bubble  splitting  noire  over 
a  range  of  velocity.  Stability  theory  is  being  used  to 
explain  the  results. 

It  can  be  concluded  that  splitting  noise  is  the  most  potent  noise 
source  when  it  occurs.  An  average  Increase  of  20  dB  over  single  phase  jet 
noise  is  common.  The  time  scale  of  each  noise  pulse  is  Independent  of 
velocity,  leading  to  the  conjecture  that  the  splitting  process  is  triggered 
by  turbulence,  but  the  resulting  unstable  motion  is  a  function  of  only 
bubble  geometry.  The  peak  sound  pressure  resulting  from  the  formation  of 
bubbles  is  dependent  on  the  details  of  how  the  bubble  was  formed.  A  simple 
model  of  this  process  correlated  well  with  the  measured  data. 

This  work  is  continuing.  Efforts  are  under  way  to  form  a  theoretical 
basis  for  the  observed  instabilities  which  lead  to  splitting.  Additional 
experiments  are  now  being  run  to  more  clearly  elucidate  the  complex  physics 
involved. 
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INTRODUCTION 


The  noise  generated  by  multiphase  flows  has  been  studied  by  many 
researchers.  Cavitation  produces  especially  high  noise  levels,  showing 
some  correlation  with  damage  resulting  from  prolonged  exposure  to  such  con¬ 
ditions.  Substantial  noise  radiation  is  also  possible  in  non-cavitating 
air-water  mixtures.  The  highest  noise  levels  result  from  the  splitting  of 
air  bubbles  when  the  two-phase  mixture  is  exposed  to  turbulent  shear  flow 
conditions. 

Hinze  [1955]  was  one  of  the  first  to  describe  the  splitting  process  of 
drops  and  bubbles  in  turbulent  flow  fields.  Many  others  have  followed, 
including  Collins  and  Knudson  [1970],  Azbel  [1981],  Sevik  and  Park  [1973] 
and  most  recently  Bentley  and  Leal  [1986]  in  a  series  of  finely  controlled 
experiments.  While  these  researchers  were  Interested  in  splitting  mecha¬ 
nics,  others  have  been  concerned  with  noise  generation  due  to  the 
splitting.  Strasberg  [1956]  was  one  of  the  first  to  postulate  on  noise 
production  due  to  many  different  fundamental  processes  in  bubble  dynamics, 
including  splitting.  Blake  [1976]  and  Gavigan  et  al.  [1974]  have  presented 
experimental  data  on  noise  generated  by  gas  jets  in  a  turbulent  wake.  In 
spite  of  the  numerous  literature  on  the  general  topic,  the  noise  associated 
with  the  splitting  of  single  bubbles  has  not  yet  been  reported.  In  addi¬ 
tion,  the  noise  due  to  formation  of  bubbles  has  only  been  given  rudimentary 
treatment.  Both  Issues  were  treated  in  this  program. 
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TEST  FACILITIES 


A  submerged  jet  facility  at  the  St.  Anthony  Falls  Hydraulic 
Laboratory,  University  of  Minnesota,  was  used  for  the  splitting  experi¬ 
ments.  A  25.4  mm  diameter  jet  issued  vertically  from  a  nozzle  into  a  tank 
of  quiescent  water.  Degassed  water  (5-7  ppm  total  dissolved  gas)  was 
pumped  through  the  closed  loop  system,  with  jet  velocities  of  up  to  20  m/s 
possible.  The  ambient  pressure  in  the  tank  could  be  controlled  between  0.1 
atm  and  1.7  atm.  The  turbulence  level  at  the  nozzle  exit  plane  was  deter¬ 
mined  to  be  0.6  percent  using  a  two  component  laser  Doppler  veloclmeter. 
Figure  1  shows  the  test  facility. 

Air  bubbles  were  introduced  into  the  flow  just  downstream  of  the  tur¬ 
bulence  management  section.  Dried,  compressed  (4-6  pslg)  air  was  passed 
through  a  27-gauge  (0.203  mm  I.D.)  hypodermic  needle,  resulting  in  bubbles 
of  1.1  ±  0.1  mm  in  diameter.  The  bubbles  were  injected  along  the  cen¬ 
terline  and  carried  with  the  flow  into  the  potential  core  of  the  jet. 
Splitting  of  the  bubbles  occurred  in  the  shear  layer  of  the  jet,  usually 
near  the  end  of  the  region  of  flow  development. 

The  acoustic  signals  were  measured  with  a  miniature  hydrophone 
positioned  in  the  tank  at  four  jet  diameters  downstream  from  the  jet  exit 
plane.  The  sound  pressure  signal  was  then  passed  through  a  low  noise 
amplifier  into  an  analog  anti-aliasing  filter  and  finally  into  a  digital 
oscilloscope.  The  oscilloscope  was  set  to  capture  the  transient  bubble 
splitting  events  at  a  data  rate  of  500  kHz.  The  waveforms  were  then 
transferred  to  a  microcomputer  for  storage  and  further  analysis. 

The  bubble  creation  experiments  were  carried  out  in  a  small  rec¬ 
tangular  tank  (508  mm  by  254  mm  by  304  mm),  filled  with  distilled  water. 
Air  bubbles  were  created  by  passing  dry,  compressed  air  through  a  submerged 
hypodermic  needle.  The  inlet  pressure  was  regulated  to  27.5  kPa  and  the 
air  flow  was  controlled  with  a  small  needle  valve.  Air  flowrates  were 
measured  with  a  rotameter,  water  temperature  and  atmospheric  pressure  were 
also  recorded.  Several  needle  sizes  were  tested:  18,  21,  25,  27,  and  30 
gauge,  corresponding  to  0.838,  0.508,  0.254,  0.203,  and  0.152  mm  I.D., 
respectively. 

The  sound  pressure  was  measured  with  a  miniature  hydrophone.  The 
amplified  output  was  passed  through  analog  anti-alalslng  filters  and  into  a 
digital  oscilloscope.  The  waveform  was  acquired  on  the  oscilloscope  and 
then  transferred  to  a  microcomputer  for  storage  and  further  analysis.  A 
schematic  of  the  experimental  setup  is  shown  in  Figure  2. 
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Figure  1 .  Schematic  of  test  apparatus 
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Figure  2.  Schematic  of  experimental  apparatus 


BUBBLE  FORMATION  NOISE 


Theory 

The  noise  generated  by  air  bubble  formation  at  a  nossle  is  charac¬ 
terised  by  a  damped  sinusoidal  oscillation.  The  sound  pressure  is  a  direct 
result  of  motion  of  the  bubble  wall.  Strasberg  [1956]  has  shown  that 
simple  volume  pulsation,  or  sero  mode  oscillation,  is  dominate  in 
generating  sound.  The  natural  frequency  of  this  mode  was  given  by  Minnaert 
[1933]. 


f  -  (3yP  /pV1/2/2*R 

O  O  D 


(1) 


where  y  -  ratio  of  specific  heats  of  the  gas  in  the  bubble, 

PQ  «  ambient  pressure, 
p  ■  density  of  the  fluid,  and 
Rb  -  bubble  radius. 

These  volume  oscillations  of  the  bubble  are  excited  by  a  change  of  pressure 
on  the  bubble,  either  internally  or  externally.  By  solving  the  linear 
equation  describing  bubble  motion, 


Me'v  +  R  v  +  Kfv(t)  -  y  )  -  0 


(2) 


with  the  initial  conditions  of  (v(t)  -  v0)  -  Av0;  and  v(0),  the  pressure  at 
a  distance  r  from  the  source  is 

poj*  -pf  6t  uf  St 

P  ■  A  e  co8(uot  -  $)  ■  P  e  cos(u>0t  -  $)  (3) 


>*/\  •  • *  v  2  1 /2 

where;  A  -  Av  f  1  +  ■  -  +  (—  ]  ,  and 

o  (o  Av  u  Av  * 

oo  oo 

A  *  dissipation  constant. 

In  Strasberg* s  model,  the  assumptions  are  that  the  bubble  grows  very 
slowly  and  the  pressure  inside  the  nozzle  is  constant  and  equal  to  Pj  - 
PQ  +  P+  ■  PQ  +  2S/Rn.  These  conditions  give  a  resulting  bubble  wall 
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velocity  of  R  «  (2/3  P+/p)l/2  .  This  value  waa  then  used  as  the  Initial 
condition  driving  the  oscillation  neglecting  Av0.  The  initial  voluae  is 
then  given  by 

*(0)  -  4*R  R2b-  4*  R2(  2/3  P+/p)  1/2  (4) 

resulting  in  a  peak  sound  pressure  of 

pp '  Vr,<*'  p„  \)1/2  <5> 

The  purpose  of  our  experiments  was  to  verify  this  theory,  in  par¬ 
ticular  the  peak  sound  pressure  level.  The  oscillation  and  damping  charac¬ 
teristics  were  also  studied. 

Results 


As  expected,  the  noise  waveforms  were  damped  sinusoidal  oscillations. 
A  typical  waveform  is  shown  in  Figure  3  and  correlated  with  bubble  motions 
observed  with  high  speed  photography.  The  time  at  which  each  photograph 
was  made  is  marked  in  letters  A  through  F  on  the  time  trace.  The  small 
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Figure  3.  Typical  sound  pressure  waveform  and  corresponding 
bubble  motion. 


cylindrical  shape  below  each  bubble  Is  the  hypodermic  needle  used  to 
produce  the  bubble.  Note  that  by  position  F  the  bubble  has  moved  a  con¬ 
siderable  distance  from  the  needle.  Also  note  that  the  peak  sound  pressure 
is  generated  while  the  bubble  is  still  spherical.  Each  transient  waveform 
was  analyzed  for  the  peak  sound  pressure  and  its  frequency  spectrum.  The 
peak  sound  pressures  generated  for  each  needle  size  at  a  distance  of  r  •  20 
mm  are  shown  in  Figure  4.  The  two  lines  indicate  an  interesting  phenomenon 
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Figure  4.  Peak  sound  pressures  for  a  range  of  needle  sices 


associated  with  the  addition  of  a  light  coating  of  oil  onto  the  needle 
surface.  A  localised  change  in  surface  tension  is  noted  by  a  slight 
Increase  in  bubble  size  for  a  given  needle  diameter,  but  perhaps  more 
interesting  is  the  substantial  drop  in  peak  sound  pressure. 

The  frequency  analysis  was  carried  out  with  an  TFT  (Fast  Fourier 
Transform)  routine  on  8192  data  points.  The  results  for  all  cases  was  a 
nearly  singular  frequency  component  at  the  natural  frequency  of  the  volume 
pulsation  of  the  bubble.  The  exponential  decay  of  the  sound  pressure  was 
best  characterised  by  the  value  of  6,  the  dissipation  constant,  given  by 
Meyer  and  Skudrsyk  [1953].  This  constant,  6  <*  0.014  +  1.1  K-05  f0,  when 

input  into  the  term  describes  the  decaying  envelops  very  well  for 
all  bubble  sizes.  Additional  tests  validated  the  1/r  dependence  on  the 
peak  sound  pressure. 

Discussion 


The  classic  theory  for  bubble  formation  at  a  nozzle,  as  given  by 
Strasberg,  does  not  give  a  fully  satisfactory  result  for  our  experimental 
data.  While  the  frequency  and  decay  match  quite  well,  the  peak  sound 
pressure  is  on  average  40  dB  below  the  value  predicted  by  Strasberg* s 
theory.  Since  frequency  and  damping  are  predicted  correctly,  the  difference 
in  the  peak  amplitude  can  be  tied  directly  to  the  initial  conditions 
Av0  and  v(0)  ,  used  in  solving  the  linear  differential  equation  for  bubble 
motion.  There  is  some  prior  evidence  that  the  value  of  v(0)  predicted 
by  hydrodynamic  theory  is  not  valid.  When  compared  to  the  measured  value 
Strasberg  presents  in  his  pape^. ,  the  predicted  value  is  one  order  of  magni¬ 
tude  higher;  however,  no  comment  is  made  on  this  discre-.  ocv. 

The  need  for  a  slightly  different  model  is  evide*  The  main  area  of 
question  is  in  the  assumption  of  constant  internal  pt  re.  The  excess 
pressure  P+  initially  has  to  be  at  least  2S/Rn;  howe  <  is  the  bubble 
grows,  the  magnitude  of  P+  does  not  necessarily  remain  con  Much  less 
pressure  is  required  for  bubble  growth  than  is  required  initiate  the 
formation.  If  P+  decreases,  so  does  the  bubble  wall  velocity  R(0),  and  in 
turn  the  peak  sound  pressure  is  also  decreased.  As  the  bubble  detaches  from 
the  needle,  the  bubble  surface  must  close  in  order  to  seek  its  minimum 
energy  state.  This  closure  imparts  some  kinetic  energy  into  the  bubble 
oscillation.  By  equating  the  work  done  to  close  the  bubble  with  the  kinetic 
energy  of  the  resulting  oscillation,  one  can  solve  for  a  new  value  of  v(0). 
of  v(0).  Figure  5  identifies  the  variable  used  in  this  formulation. 

The  work  done  in  closing  the  bubble  can  be  given  as 


W  -  M9  -  2wR  (P.  -  P  )  f”  8  ds  0  •  »R3  P*0  (6) 

n  l  o  n 

o 


where  Pi  -  P0  •  P*  ■  2S/Ri,.  The  kinetic  energy  of  the  oscillation  is  given 
by: 


Ek  -  1/20^  +  mad<J)R2 


(7) 
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By  equating  equations  6  and  7,  and  solving  for  the  volume  velocity,  we  get 
3/2  1/2 

v(0)  -  4irRn  (S0/p)  .  Using  this  Initial  condition,  new  values  of  peak 

sound  pressure  can  be  calculated,  which  are  of  the  same  order  as  our 
measurements.  The  results  are  summarized  below,  again  for  r  ■  20  mm  and 
without  oil. 


Needle  size 

(mm  I.D.) 

Bubble  size 

(mm  dla.) 

P 

meas. 

Pmodel 
(dB  re  lyPa) 

p 

theory 

0.838 

3.812 

mm 

153.5 

179.4 

0.508 

3.048 

148.9 

179.6 

0.254 

2.370 

142.1 

180.5 

0.203 

2.052 

140.4 

180.2 

0.152 

1.778 

139.1 

137.9 

180.2 

j 

j 
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SPLITTING  NOISE 


Experiments 

The  splitting  noise  experiments  were  Initially  designed  utilizing 
critical  Weber  number  criteria  developed  by  Hlnze  [1955]  and  further  tested 
for  a  submerged  jet  by  Sevik  and  Park  [1973] .  For  a  given  bubble  Injector 
diameter,  It  is  possible  to  go  from  a  non-splitting  to  a  splitting  regime 
by  simply  varying  the  jet  velocity.  Figure  6  shows  this  concept  for  a  0.1 
mm  diameter  bubble  Injector  and  a  25.4  mm  diameter  jet.  The  details  of  the 
calculation  can  be  found  in  Frizell  [1987].  The  data  presented  in  this 
paper  cover  jet  velocities  from  2.7  m/s  to  17.5  m/s.  Bubble  splitting  was 
observed  at  all  velocities  tested.  Single  phase  jet  noise  levels  along 
with  fifteen  splitting  events  were  recorded  at  each  velocity.  Photographs 
were  taken  to  document  each  test  run.  The  results  are  discussed  below. 

Discussion 

The  sound  pressure  waveform  of  a  splitting  event  is  characterized  by  a 
dumped,  multi-frequency  oscillation  (see  Figure  7).  In  most  cases,  the 
results  Indicated  that  there  was  some  type  of  volume  oscillation  present, 
although  the  frequencies  of  the  oscillations  generally  predicted  a  bubble 
size  closer  to  the  initial  bubble  diameter  rather  than  the  size  resulting 
after  splitting.  A  comparison  of  the  acoustic  intensity  in  1/3  octave 
bands  shows  a  general  increase  in  level  and  a  broading  of  the  spectra  at 
the  higher  frequencies  (see  Figure  8). 

Some  of  the  broadening  effects  can  be  explained  by  the  measurement 
position.  At  5  kHz,  for  example,  the  hydrophone  is  positioned  in  the 
acoustic  near  field.  On  the  other  hand,  the  measurement  position  is 
marginally  in  the  acoustic  far  field  at  50  kHz.  The  acoustic  intensity  is 
expected  to  vary  with  u  in  the  near  field  and  to  a  higher  power  of  U  in 
the  acoustic  far  field,  hence  a  possible  reason  for  the  broadening  of  the 
spectrum  with  increasing  velocity.  Since  the  bulk  of  the  acoustic  energy 
is  centered  about  frequencies  corresponding  to  the  acoustic  near  field,  the 
p^ak  pressure  follows  a  TJ2  relationship  (U**  for  intensity)  as  compared  to 
IT  for  Blakes  [1976]  measurements  that  were  made  in  the  acoustic  far  field. 
This  dependence,  along  with  a  comparison  to  single  phase  jet  noise  and 
bubble  creation  noise,  is  shown  in  Figure  9. 

Further  analysis  of  the  splitting  waveforms  yields  interesting 
information  about  the  characteristic  time  scale  of  the  events.  If  the 
splitting  is  to  be  an  event  dominated  by  the  turbulent  field,  one  would 
expect  that  the  time  scale  of  the  splitting  would  be  of  the  same  order  as 
the  characteristic  time  scale  of  the  turbulence.  However,  it  was  found 
that  the  time  scale  of  the  splitting  was  three  orders  of  magnitude  smaller 
than  that  of  the  turbulence.  This  comparison  is  made  using  time  scales 
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j  Figure  6.  The  solid  line  in  the  figure  is  critical 

Weber  number  versus  velocity.  The  dotted 

!  line  is  a  prediction  of  bubble  size  for 

i  an  0.1  mm  injector. 

I 
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Figure  9.  Comparison  of  peak  sound  pressure  in  a  jet 
with  and  without  bubbles . 
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typical  of  turbulent  fluctuations  that  are  spatially  coherent  over  one 
bubble  diameter.  Following  Sevik  and  Park  [1973],  the  following  definition 
is  used: 


.2  1/3 

rf-> 


(9) 


where  e  Is  the  dissipation  rate  and  d  is  the  bubble  diameter.  In 
addition,  the  splitting  time  scale  did  not  vary  with  jet  velocity  (see 
Figure  10).  This  fact  may  lead  to  a  new  postulation  that  while  the 
turbulence  is  responsible  for  exciting  the  bubbles  into  oscillations, 
resulting  in  severe  deformation,  the  actual  splitting  is  the  direct  result 
of  an  instability. 
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CONCLUSIONS 


The  results  of  the  splitting  experiments  have  shown  that  the  noise 
generated  by  single  bubbles  splitting  in  a  turbulent  shear  flow  can  be 
substantial.  An  average  rise  of  20  dB  above  single  phase  jet  noise  is 
common. 

Comparison  of  the  time  scale  of  a  splitting  event  with  typical 
turbulence  time  scales  leads  to  the  conjecture  that  splitting  is  the  result 
of  an  instability  that  is  triggered  by  turbulence,  but  the  final  stages  of 
the  process  are  a  function  of  bubble  geometry. 

The  formation  of  a  bubble  at  a  nozzle  generates  a  sound  pulse  in  the 
form  of  a  damped  sinusoidal  oscillation.  The  frequency  of  oscillation  is 
associated  with  the  volume  pulsation  of  the  bubble.  The  damping 
characteristics  in  the  frequency  range  studied  are  dominated  by  acoustic 
and  thermal  damping,  and  are  well  defined  by  theory.  The  peak  sound 
pressure  is  a  function  of  the  method  used  to  form  the  bubble,  the  initial 
wall  velocity  being  the  most  important  parameter.  The  values  predicted  by 
Stranberg's  model  did  not  match  our  experimental  data.  The  assumptions 
were  reviewed  and  a  simple  model  formulated  to  predict  the  initial  volume 
velocity.  The  results  of  this  model  predicted  our  experimental  data 
usually  within  3  to  6  dB.  The  Influence  of  a  coating  of  oil  on  the  needle 
surface  is  noticed  in  two  ways:  1)  a  change  in  the  localized  surface 
tension,  resulting  in  a  slightly  larger  bubble  than  for  the  nonollled 
needle,  and  2)  approximately  a  10  dB  drop  in  the  peak  sound  pressure  level. 
This  drop  in  noise  has  not  been  explained  as  of  yet,  but  opens  an  area  for 
further  research. 
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